Derivation of Maxwell-Bloch-type equations by projection of quantum models 
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A simple algebraic procedure is described for deriving Maxwell-Bloch-type equations from single- 
atom cavity quantum electrodynamics (cavity QED) master equations via orthogonal projection onto 
a manifold of semiclassical states. In particular the usual Maxwell-Bloch Equations are obtained — 
up to a state-dependent correction factor of order unity — straightforwardly from the unconditional 
Jaynes-Cummings master equation. The technique of projecting onto a semiclassical manifold can 
also be applied with conditional master equations (quantum filters) , leading to stochastic simulation 
models that include multiplicative noise terms associated with fluctuations of the atomic dipole. The 
utility of such models is briefly explored in the context of single-atom absorptive bistability. 

PACS numbers: 42.50.Pq,42.50.Lc,42.65.Pc,03.65.Sq 



For many years, cavity quantum electrodynamics (cav- 
ity QED) in the optical regime has provided a canon- 
ical setting for theoretical and experimental investiga- 
tions of non-equilibrium quantum statistical mechanics 
and the quantum-classical transition [3, 0, Histori- 
cally much attention has been focused on the comparison 
of semi-classical models of cavity QED, exemplified by 
the Maxwell-Bloch Equations (MBE's) and quantum- 
mechanical models such as the Jaynes-Cummings master 
equation It is now generally accepted that the MBE's 
are most directly applicable in scenarios where many in- 
tracavity atoms each contribute weakly to the overall co- 
operativity [H[E,@, whereas quantum- mechanical models 
are required to account for experimental data in single- 
atom scenarios with strong coupling 0, Q. There are 
however some subtleties to this distinction as, e.g., it can 
be shown that MBE's govern the dynamics of mean val- 
ues of quantum observables in single-atom cavity QED 
scenarios in the "bad cavity limit" with weak excita- 
tion 0. Rigorous distinction between semi-classical and 
quantum phenomena in cavity QED is of fundamental 
interest for both the study of mesoscopic physics and 
quantum information science 

In investigations based on numerical simulation of 
quantum models [IE GH d GS d, 

it has been ob- 
served that the MBE's actually provide fairly accurate 
guidance when searching parameter space for bifurcation- 
like phenomena in single-atom cavity QED with strong 
driving fields. This extended relevance of the MBE's 
in regimes of strong coupling and nonlinear dynamics is 
somewhat surprising, and there has been no prior anal- 
ysis of the quantitative relation between the MBE's and 
any fully quantum-mechanical model of single-atom cav- 
ity QED that would be valid under such conditions. In 
this article I point out that a simple technique, intro- 
duced recently into quantum optics by Van Handel [HI], 
can be used to show that the MBE's essentially corre- 
spond (up to a factor of order unity) to an orthogonal 



projection of the single-atom cavity QED master equa- 
tion onto an intuitive manifold of semiclassical states. 
No assumptions regarding parameter regime are required. 
Applyin g th e same technique to stochastic master equa- 
tions 

H HE E3, HI (quantum filters for detection 
of the atomic fluorescence [IE HH, I arrive at modified 
MBE's with noise terms that constitute semi-classical 
cavity QED models incorporating atomic dipole fluctua- 
tions. These models provide an interesting complement 
to stochastic extensions of the MBE's that have been pro- 
posed for incorporating electromagnetic field fluctuations 

As usual we consider cavity QED models in which a 
two-level atom is coupled to a single mode of the elec- 
tromagnetic field [J, UK]. Our general strategy [lj| for 
obtaining semi-classical MBE's is to project a quantum 
equation of motion, which we denote abstractly 



dd t = C [0 t 



(1) 



onto the manifold of states 
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where \a) is a coherent state with complex amplitude 
a for the electromagnetic field mode, j± is the atomic 
dipole decay rate and go is the "vacuum Rabi frequency" 
characterizing the strength of the atom-cavity coupling. 
In our parametrization <J x ^ VtZ are the usual Pauli opera- 
tors and / represents the identity operator on the atomic 
Hilbert space. We are thus restricting the system dynam- 
ics to unentangled quantum states in which the field is 
semi-classical and the atom assumes an arbitrary mixed 
state. Our chosen parametrization of the atom-cavity 
density operator in terms of the real scalar variables p r , 
Pi, D, x r and Xi will lead us to Maxwell-Bloch-type equa- 
tions. The extra parameter n is included for technical 
reasons to allow us to enforce normalization of the den- 
sity operator in the projected equations, and will drop 
out of the final equations. 
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Orthogonal projection of Eq. JT]) is accomplished via 
the general rule 



<h>< ( ^-,d9t^ ' d Pi = \§ P ~ ,< ^ t 



(11) 



ll S pan{ ! / i }Fj = A 



(3) 



and 



with x «-> <i#t and with {^} corresponding to a full set of 
tangent vectors to the parametrized manifold of states: 

s v .\ <_> {®R ®£. ®E ®£. J^R. dp \ (4) 

1 ' 1 \ dn ' dp r ' ' dD ' 9a; r ' do^ / 

We choose the inner product (X, Y) = Tr[X*Y] on the 
space of atom-cavity density operators, which will prove 
to be especially convenient because Pauli operators play 
a central role in our parametrization. Note that other 
choices of inner product could be made and might lead 
to somewhat different final equations. 

After some straightforward calculations we find 

dp 1 T i \ / i dp 1 I i / I , rS 
— = -I®\a)(a\, —=--<T g ®\a){a\, (5) 



dp 1 



dp 1 



— = --^a,9\a)(a\, ^ = ^®WR (6) 
and 

Q 

-7^- = Pat ® [(a* - a*)\a)(a\ + \a)(a\(a- a)] , 



dp_ 

dxi 

where 



Pat ® [i(a* - a*)\a)(a\ - i\a)(a\(a - a)] , (7) 



= i (nl - V2p r (T x + V2pi<r y - Ddz) . (8) 



Pat 



As given these tangent vectors are mutually orthogonal 
but the last two arc not normalized: 



dp dp \ / dp dp 

dx r ' dx r / \ dxi ' dxi 



= n 2 + 2p 2 r + 2p 2 + D 2 . (9) 

This provides the denominator required in Eq. 
Finally we simply write 

dp t = H spaa { Vi }[ddt\ 

d P j , d P a , d P A 
= —dn + — dp r + ^—dp l 

dn dp r dpi 

■ dp j „ dp dp 
+ ^-=;dD + - — dx r + ^—dxi, 
dD dx r dxi 

and thus, using Eq. © and orthogonality of the tangent 
vectors, 



dp 



dn=(^-,dO t ), dD = (^,d8 t 



l dp 

= \dD^' 



(10) 



dx r = {n 2 + 2p 2 r + 2p 2 + D 2 )- 1 ( 4?-,d8 t 

\dx r 

,2 , r>2 I n„2 , r>2\-l / dP 



d Xi = (n 2 + 2p 2 + 2pf + D 2 )-^^,de t J. (12) 

All that remains in order to obtain concrete semi-classical 
equations of motion for the real scalar parameters is to 
insert a specific quantum model for £[#(]. It is interesting 
to note that, in any given time-step, the error associated 
with projection onto the target manifold can be quanti- 
fied by computing the norm of (d0 t — dpt). 

We first consider the unconditional Jaynes-Cummings 
master equation [H, written in a rotating frame at the 
frequency lui of the external driving field (h = 1): 

d0 t = -i[H, 6 t )dt + n(2a6 t a* - a*a6 t - 9 t a*a)dt 

+7 ± (2cr(9 t cr* - a*aB t ~ t a*a)dt, 
H = A c a* a + A a a* a + igo(a* a — aa*) 

+i£(a*-a). (13) 

Here k is the cavity field decay rate, a is the field annihila- 
tion operator, a is the atomic dipole (lowering) operator, 
A c = uj c — uji and A a = uj a — uji are detunings of the cavity 
and atomic resonance frequencies, and £ represents the 
complex amplitude of the coherent driving field. With a 
view towards matching the usual convention for MBE's, 
we apply the scalings 

P«-— 7=> D*--D, «<-7±i, (14) 



where n = 1±_ /^9o ^ s ^ ne critical photon number. 

After some rather tedious calculations, and adopting 
the conventional parameter definitions A = A a /j±, = 
A c /k, k = k/j±, C = 5q/2k7^, y — E/k^/tiq, we arrive 
at the dimensionless equations 



dn dD 
~dt ~ ' ~dt 



-2(D-l)-2(p r x r +p iXi ), (15) 



= -pr + Api + DXr, = -pi-Ap r + DXi, (16) 

dt dt 



dt 

dxi 

~dt 



= -k {x r - Qxi - Re[y] + 2Cp r F} , 
= -k {xi - <dx r - Im[y] + 2CpiF} , 



F = 



l + 2p 2 + 2p 2 + D 2 • 



(17) 



We have invoked the fact that n is constant and set n = 1 
in the equations above (as is required for normalization 
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of the projected density operator p t ). Looking at the re- 
sulting equations, we see that we would exactly recover 
the usual MBE's by setting F — > 1. Given the definitions 
and scalings used above, we can write F — 2/(1 + IS^ 2 ), 
where S is the atomic Bloch vector. It follows that for 
atomic pure states F — ► 1 and we can justify the usual 
MBE's if we have reason to believe that this purity is 
nearly preserved by the dynamics (e.g., under very weak 
excitation). Likewise we see that I < F < 2, but compar- 
ison of Eqs. (fl"5]) - (fTT)) with the usual MBE's in a standard 
setting such as absorptive bistability shows that even this 
factor can lead to significant corrections (see below). 

Turning now to the projection of conditional master 
equations, we first note that terms associated with condi- 
tioning upon homodyne/heterodyne detection of the cav- 
ity output field vanish upon projection. We therefore be- 
gin with the quantum filter for homodyne detection of the 
atomic fluorescence (2~j| . written in linear Stratonovich 
form as appropriate for the projection method (as dis- 
cussed in [la]): 



dO t = -i[H, 6 t ]dt + K(2ad t a* - a*a6 t 
-j ± ((T*(r8 t + 9 t a*a)dt 
+ y/2jl{a0t + 9 t (7*)ody t . 



t a*a)dt 



(18) 



Here dyt represents the measured photocurrent that 
drives the filter. Here we are not yet scaling time by 
7x- The parameters and the Hamiltonian are as in the 
unconditional master equation (113|) . 

Projection leads to equations including a nontrivial 
evolution of the normalization, dn — j±(D — n)dt — 
2^/7j_p. r o dy t , which is not surprising since we started 
from an unnormalized model. We therefore transform to 
normalized variables by defining p r = p r /n, pi = Pi/n, 
D = D/n. Converting the resulting equations to Ito form 
we obtain the projected filter 



dp r 
dpi 
dD 



dx, 
If 
dx 
~dt 



l±{-ipr + Ap t + 2p r D + Dx r + Apl)dt 
+V7l(2p 2 + D- l)dy t , 
7_i_(-Ap r - pi + Dxi + 4plpi)dt 

+2y/~fIPrPidyt, 

7 ± (2 - 2D - 2p r x r - 2piXi - 4p 2 . + Ap 2 r D)dt 
-2Vl±(pr -p r b)dy t , 

—k^±(x r — Qxi — Re[y] + 2Cp r F), 



—1 = -k-y±(xi + @x r -Im[y]+2CpiF). 



(19) 



As it is our aim to obtain stochastic equations that 
could be used for Monte Carlo simulation of semi-classical 
cavity QED dynamics, we now make the substitution 



dy t — > dW t - 2^/jTp r dt, 



(20) 



where dWt is an Ito increment and —2y/j±p r represents 
an approximation within the semi-classical state space 
of ^2j± Tr [a x 6t] /Tr [9t] , which is the expected value of 



the measured signal dyt- Making this substitution, and 
scaling time by t <— j±t, dW t <— y/j±dWt, we obtain 



dp r = (-p r + Apt + Dx r )dt + (2p 2 r + D - l)dW t , 

dp t = (— pi — Ap r + Dxi)dt + 2p r pidWt, 

dD = (2 -2D -2p r x r -2piXi)dt + 2p r (D - I) dW t , 

Re[y] + 2Cp r F), 
®x r - Jm[y] + 2CpiF). 
Using these equations it can be shown that 



— — — k(x r Qxi 
dt 

dxi , 



(21) 



d(2p 2 r 



2p 2 l+ b 2 )tx(2pl + 2pl 



D - 1) 



(22) 



and thus vanishes under the initial condition 2p 2 + 2p 2 + 
D 2 = 1. Hence for an initial atomic pure state we can 
simplify the field evolution equations to 



dt 
dxi 
~dt 



-k(x r — Qxi — y + 2Cp r ), 
-k(xi + Qx r + 2Cpi). 



(23) 



Following an analogous procedure for a linear 
Stratonovich quantum filter for heterodyne detection of 
the atomic fluorescence, 

d9 t = -i[H, 9 t ]dt + n(2a9 t a* - a*aO t - 9 t a*a)dt 
-j±(<T*a9t + 9 t a*a)dt 
+ ^I(a9t + 9 t a*)oRe[dyt} 
+i^2^Z(<79t~9 t a*)olm[dy t ], (24) 

where dyt now represents a complex photocurrent as in 
[H,[i3|, we obtain 

dpr = (—pr + Api + Dx r )dt 

+ -^(2p 2 r + D- l)dW[ - V2p r pidWi, 
v2 



dpi 



dD = 



dt 

dxi 

~dt 



(—pt — Ap r + Dxi)dt 

+V2p r pidW[ - -j={2p: 
v2 



-2 



D - l)dWl 

(2 — 2D — 2p r x r — 2piXi)dt 

+V2p r (b - l)dW[ - V2pi(D - l)dW l t , 

-k(x r - ®Xi - Re[y] + 2Cp r F), 



= -k(xi + Qx r - lm[y] + 2Cp l F). 



(25) 



Here dW[ and dW\ are independent Ito increments. 

Summarizing the analytic results, we have shown that 
projection according Eq. ([3]) of quantum master equa- 
tions onto the semi-classical manifold of states speci- 
fied by Eq. ([2]) leads to Maxwell-Bloch-type equations 
for single- atom cavity QED. Starting from the uncondi- 
tional master equation fl~3")) we obtain the model (jTSJ)- 
(I17p . From the quantum filter for homodyne detection of 
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FIG. 1: (a) Simulated trajectory of x r from model (|25|l with 
parameters of absorptive bistability. (b) Joint distribution 
of x r and Xi from the same simulation, (c) Comparison of 
equilibrium solutions according to the usual MBE's (dashed) 
and our projected equations (|15p ~ (|17p . 

the atomic fluorescence, Eq. (fl~8|). we obtain the projected 
filter (|19[) and then use the approximation (f2"0")) to obtain 
the stochastic simulation model (f2"Tj) . From the hetero- 
dyne quantum filter, Eq. (|24[) . we obtain by an analogous 
procedure the stochastic simulation model (f25|) . Condi- 
tioning terms associated with homo/heterodyne detec- 
tion of the cavity output field disappear upon projection. 

We turn finally to a brief examination of the behav- 
ior of our heterodyne model (|25p in a parameter regime 
for which the usual MBE's exhibit absorptive bistabil- 
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